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Running WRF-Sfire with real data^ 

Process of running real WRF-Sfire case is very similar to running the 
regular WRF and consists of the following steps: 

1. Getting needed data 

2. Defining the domain setup in namelist . wps 

3. linking meteorological files in grib format using linkgrib . sh 

4. Uncompressing (ungribing) meteorological files using 

ungrib . exe 

5. Preparing high-resolution fire data using convert_geotif f 

6. Preparing geo_em.dOX files containing all static data for the 
WRF domains using geogrid.exe 

7. Interpolating meteorological data to WRF domains defined by 

geo_em . dOX files using metgr id . exe 

8. Defining simulation setup in namelist . input and 
namelist . fire and creating wrfinput_dOX and 
wrfbdy dOX filesusingreal.exe 




Running WRF-Sfire with rea 
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Obtaining data for geogrid 


• Geogrid needs the static data available from: 
http://www2.mmm.ucar.edu/wrf/src/wps files/geog v3.1.tar.gz 

• This dataset must be extracted (tar -xvf geog_v3.1.tar.gz) 

• Its location must be specified as geog_data_path in 
namelist . wps 

for instance: 

geog_data_path = /home/ j smith/ geog/ 

• additionally high resolution elevation data and fuel data are needed 
in geotiff format 

• for US they are available from: http://landfire.cr.usgs.gov/viewer/ 



Obtaining fire data for geogrid 



http://landfire.cr.usgs.gov/viewer/ 
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Obtaining fire data for geogrid 


http://landfire.cr.usgs.gov/viewer/ 
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Obtaining fire data for geogrid 
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http://landfire.cr.usgs.gov/viewer/ 








Obtaining fire data for geogrid 
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l3 us_130 13 Fire Behavior Fuel Models-Anderson 

GeoTIFF_with_attribs i ] 
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□ us_130 40 Fire Behavior Fuel Mod els -Scott/ Burg an 

Not selected. 

□ us_130 Biophysical Settings 

Not selected. 

□ us_130 Environmental Site Potential 

Not selected. 

□ us_130 Existing Vegetation Cover 

Not selected. 

□ us_130 Existing Vegetation Height 

Not selected. 

□ us_130 Existing Vegetation Type 

Not selected. 

□ us_130 Fire Regime Groups 

Not selected. 

□ us_130 Forest Canopy Base Height 

Not selected. 

□ us_130 Forest Canopy Bulk Density 

Not selected. 

□ us_130 Forest Canopy Cover 

Not selected. 

□ us_130 Forest Canopy Height 

Not selected. 

□ us_130 LCP 13 Fire Behavior Fuel Models-Anderson 

Not selected. 

□ us_130 LCP 40 Fire Behavior Fuel Models-Scott/Burgan 

Not selected. 

□ us_130 Mean Fire Return Interval 

Not selected. 

□ us_130 Percent Low-severity Fire 

Not selected. 

□ us_130 Percent Mixed -severity Fire 

Not selected. 

□ us_130 Percent Replacement-severity Fire 

Not selected. 

□ us Aspect 

Not selected. 

$ us Elevation 

| CeoTIFF : ] 

IZIP ;] 

| HTML * ) 

□ us Slope 

Not selected. 

□ Fire Planning Units 

Not selected. 

□ LANDFIRE Zones 

Not selected. 

Q LANDFIRE ECOMAP Subsections Map 

Not selected. 


|| Delivery Options: 

||Maximum size (MB) per piece: | 250 i] 





Converting fire data to WPS form 

• Conversion of the fuel data and high resolution topography can 
be performed using convert_geotiff program available from 
https://github.com/ibeezley/convert geotiff 

• Conversion of the fuel data provided in geotiff format can be 
performred using the following command: 

. /convert_geotif f . x -c 13 -w 1 -u "fuel category" -d 

"Anderson 13 fire behavior categories" ../your fuel. tiff 


• Conversion of the high resolution topogrpahy provided as a 
geotiff file can performes similarily using the following 
command: 

. /convert_geotif f . x -u meters -d 'National Elevation Dataset 

1/3 arcsecond resolution' . . /your_elevation . tif f 

• note that convert_geotif f . x should be executed from the 
locations where the fuel and elevation data are to be stored 


flp 

Converting fire data to WPS formal 

• convert_geotiff works should support all WRF projections but 
lat-lon is the safest choice EPSG 4326 WGS84. 

• Convert_geotiff will generate an index file and a set of files in 
WPS intermediate format containing the fuel and elevation data 
ready to be processed by geogrid.exe 

• Here is an example of the content of the folder where 
convert_geotiff generated a set of WPS intermediate files: 



Converting fire data to WPS formal^ 

The index file for fuel should look like below: 

projection = albers_nad83 

truelatl = 29.500000 

truelat2 = 45.500000 

stdlon = -96.000000 

known_x = 1 

known_y = 606 

known_lat = 39.747818 

known_lon = -107.373398 

dx = 3.000000e+01 

dy = 3.000000e+01 

type = categorical 

signed = yes 

units = "fuel category" 

description = "Anderson 13 fire behavior categories" 

wordsize = 1 

tile_x = 100 

tile_y = 100 

tile_z = 1 

category_min = 1 

category_max = 14 

tile_bdr = 3 

missing_value = 0.000000 
scale_factor = 1.000000 
row_order = bottom_top 
endian = little 



The physical domain is configured in the geogrid section of namelist . wps in 
the WPS directory. In this section, you should define the geographic projection 
with map_pro j , truelatl, truelat2, and stand_lon. Available 
projections include ' lambert ' , 'polar', ' mercator', and 'lat-lon'. 


The lower left corner of the domain is located at ref Ion longitude and 
ref lat latitude. The computational grid is defined~by e_we/e_sn, the 
number of (staggered) grid points in the west-east/south-north direction, and the 
grid resolution is defined by dx and dy in meters. 


We also specify a path to where we will put the static dataset that geogrid will 
read from, and we specify the resolution of the data we plan to use. 


Sgeogrid 

e_we 

e_sn 

geog_da ta_r e s 
dx = 60 , 


43, 

43, 

'30s' , 


dy = 60 , 


map_pro j = 
ref_lat = 
ref_lon = 
truelatl = 
truelat2 = 
stand Ion = 


' lambert' , 
39.70537, 
-107.2907, 
39.338, 
39.338, 
-106.807, 


geog_data_j?ath = /wrf data/geog ' 

/ 


configuration in namelist . input 

&time control 


run days 

= 

o. 



run hours 

= 

0, 



run minutes 

= 

20, 



run seconds 

= 

0, 



start year 

= 

2006, 

0001, 

0001, 

start month 

= 

02, 

\ — 1 
o 

01, 

start day 

= 

23, 

< 

\ — 1 
o 

01, 

start hour 

= 

CM 
\ — 1 

o 

o 

o 

o 

start minute 

= 

43, 

< 

\ — 1 
o 

01, 

start second 

= 

00, 

< 

o 

o 

< 

o 

o 

end year 

= 

2006, 

0001, 

0001, 

end month 

= 

02, 

01, 

01, 

end day 

= 

23, 

< 

\ — 1 
o 

01, 

end hour 

= 

13, 

< 

o 

o 

< 

o 

o 

end minute 

= 

O 

o 

600, 

600, 

end second 

= 

0, 

< 

o 

o 

< 

o 

o 

history interval s 

= 

5, 

30, 

30, 

frames per outfile 

= 

1000, 

1000, 

1000, 

restart 

= 

. false . , 


restart interval 

= 

5 



io form history 

= 

2 

▲ 

▲ 

io form restart 

= 

2 

t 

t 

io form input 

= 

2 

1 

1 

io form boundary 

= 

2 



debug level 

= 

i d01 d02 

d03 





configuration in namelist . input 


&doma±ns 





time step 

= 

0, 



time step tract num 

= 

3, 



time step tract den 

= 

10, 



max dom 

= 

1, 



s we 

= 

1, 

1, 

1, 

e we 

= 

20, 

43, 

43 

s sn 

= 

1, 

1, 

1, 

e sn 

= 

32, 

43, 

43 

s vert 

= 

1, 

1, 

1, 

e vert 

= 

41, 

41, 

41 

dx 

= 

o 

LO 

30, 

10 

dy 

= 

50, 

30, 

10 

ztop 

= 

600, 

1500, 

1500, 

grid id 

= 

1 , 

2, 

3, 

parent id 

= 

0, 

1, 

2, 

i parent start 

= 

0, 

1, 

1, 

j parent start 

= 

0, 

1, 

1, 

parent grid ratio 

= 

1, 

2, 

3, 

parent time step ratio 

= 

1, 

2, 

3, 

feedback 

= 

1 , 



smooth option 

= 

0 



sr x 

= 

10, 

o. 

0 

sr y 

= 

10, 

o. 

0 



configuration in namelist . 


input 


fiphysics 

mp_physics 

ra_lw_physics 

ra_sw_physics 

radt 

sf_sfclay_physics 

sf_surface_j?hysics 

bl_pbl_physics 

bldt 

cu_physics 

cudt 

isf fix 

if snow 
icloud 

num_soil_layers 

mp_zero_out 

/ 


= 1 , 1 , 1 , 

= 1 , 1 , 1 , 

= 1 , 1 , 1 , 

= 30 , 30 , 30 , 

= 1 , 1 , 1 , 

= 1 , 1 , 1 , 

= 1 , 1 , 1 , 

= 0 , 0 , 0 , 

= i, i, 0 , 

= o, 0, 0, 

= 1 , 

= 1 , 

= 1 , 

= 5 , 

= 0 , 


/ 


configuration in namelist . input 


fidynamics 


rk ord 

= 

3, 



diff opt 

= 

2, 



km opt 

= 

2, 



damp opt 

= 

2 , 



zdamp 

= 

5000. , 

5000. , 

5000. , 

dampcoef 

= 

0.2, 

0.2, 

0.2 

khdif 

= 

< 

LO 

O 

O 

LO 

O 

O 

LO 

O 

o 

kvdif 

= 

< 

LO 

O 

• 

o 

< 

LO 

O 

• 

o 

0 . 05, 

smdiv 

= 

0 .1, 

0 .1, 

0 .1, 

emdiv 

= 

< 

\ — 1 
o 

o 

< 

\ — 1 
o 

o 

o 

o 

> 

epssm 

= 

< 

\ — 1 

• 

o 

< 

\ — 1 

• 

o 

0 . 1 

mix full fields 

= 

. true . , 

. true . , 

. true . , 

non hydrostatic 

= 

. true . , 

. true . , 

. true . , 

h mom adv order 

= 

5, 

5, 

5, 

v mom adv order 

= 

3, 

3, 

3 , 

h sea adv order 

= 

5, 

5, 

5 ' Tracer for smoke 

v sea adv order 

= 

3, 

3, 

3 , 

time step sound 

= 

20, 

20, 

20 , representation 

moist adv opt 

= 

1, 

1, 

1 f (requires 

scalar adv opt 

tracer opt 

= 

1, 

2, 

1, 

2, 

0 ' t jiamelist.fire_emissi 




configuration in namelist . input 


&bdy_control 

periodic x 

= 

. false . , 

. false . , . false . 

r 


symmetric xs 

= 

. false . , 

. false . , . false . 

r 


symmetric xe 

= 

. false . , 

. false . , . false . 

r 

settings 

open xs 

= 

. false . , 

. false . , . false . 

r 

open xe 

= 

. false . , 

. false . , . false . 

r 

for 

periodic y 
symmetric ys 

= 

. false . , 
. false . , 

. false . , . false . 
. false . , . false . 

r 

r 

ideal 

symmetric ye 

= 

. false . , 

. false . , . false . 

r 

cases 

open ys 

= 

. false . , 

. false . , . false . 

r 


open ye 

= 

. false . , 

. false . , . false . 

r 


spec bdy width 

= 

5 , 




spec_zone 

= 

1 , 




relax zone 

= 

4 , 




specified 

= 

. true . , . 

false . , . false . , 



nested 

= 

. false . , 

. true . , . true . , 




/ 

/ 



configuration in namelist . input 


&fire 

ifire 

fire_fuel_read 

f ire_fuel_cat 

! ignition 



Fire model will be run in d03 in this case 


! integer, 
! integer, 
! integer. 


= 0: no fire, = 2: SFIRE 

-1: read in from wrf input 

if specified which fuel category? 


fire_num_ignitions = 
fire_ignition_start_lonl= 
fire_ignition_start_latl 
fire_ignition_end_lonl = 
fire_ignition_end_latl = 
fire_ignition_radiusl = 
f i re_igni tion_s tar t_time 1 
f i re_igni tion_end_time 1 

t 

I Ignition parameters set for the fire domain only 


o, o, 1, 

0 , 0 ,- 107 . 293 , 
0 , 0 , 39 . 6986 , 
0 , 0 ,- 107 . 293 , 
0 , 0 , 39 . 7109 , 
0 , 0 , 18 , ! 

0 , 0 , 2 , ! 

0 , 0 , 2 , 


integer, only the first f ire_num_ignition used, up to 5 
start points of ignition lines, longitude 
start points of ignition lines , latitude 
end points of ignition lines 
end points of ignition lines 
all within this radius (m) will ignite, > fire mesh step 
sec for ignition from the start 


sec for ignition from the start 


fire_topo_from_atm= 0, 0, 0, ! 0 = expect fire mesh topo , 1 = from atmosphere 

fmoist_run = . false . , . false . , . true . , 

fmoist_interp . false . , . false . , . true . , 

f ire_fmc_read = 0, 0, 0, ! 0 not set use wrf input, 1 from namelist . fire, 2 read from file 

in ideal ! 



Setting up geogrid to process fire data 


All the static data being processed by geogrid.exe are 

defined in GEOGRID . TBL 

In order to enable fire data processing in geogrid.exe : 

1 . GEOGRID . TBL should be linked to GEOGROD . TBL . FIRE 
In -s GEOGRID. TBL. FIRE GEOGRID. TBL or copied 
to WPS directory 

2. Fuel and high-resolution height section of GEOGRID . TBL 
have to be edited so the paths correspond to the actual 
locations of fuel and elevation data processed using 
convert geotiff (in our case ned_data, landfire_data) 


m 

Setting up geogrid to process fire datc^ 

Editing GEOGRID. TBL (WPS/ geogrid/ GEOGRID . TBL) 



name=NFUEL_CAT 

priority=l 

dest_type=categorical 
dominant_only=NFUEL_CAT 
z_dim_name=f uel_cat 
halt_on_missing=yes 

interp_option=def ault : nearest_neighbor 

abs_path= . /path_to_your_f uel_data 

subgrid=yes 


name=ZSF 

priority = 1 
dest_type = continuous 
df_dx=DZDXF 
df_dy=DZDYF 

smooth_option = smth-desmth_special ; smooth_passes=l 
halt_on_missing=yes 

interp_option = default : average_gcell ( 4 . 0 ) +four_pt+average_4pt 

abs_path= . /path_to_your_elevation_data 

subgrid=yes 


Domain configuration in WPS 

The share section of the namelist . wps defines the fire subgrid refinement in 
subgrid_ratio_x and subgrid_ratio_y. This means that the fire grid 
will be a 10 time refined grid at a resolution of 6 meters by 6 meters. The 
start_date and end_date parameters specify the time window that the 
simulation will be run in. Atmospheric data must be available at both temporal 
boundaries. The interval_seconds defines the time interval between 
atmospheric data. For our example, we will be using the NARR dataset which is 
released daily every three hours or 10,800 seconds. 

&share 

wrf_core = 'ARW' , 
max_dom = 1 , 

start_date = ' 2005-08-28_12 : 00 : 00 ' , 
end_date = ' 2005-08-28_15 : 00 : 00 ' , 
interval_seconds = 10800, 
io_f orm_geogrid = 2 , 
subgrid_ratio_x = 10, 
subgrid_ratio_y = 10, 

/ 



£>jr rJt 


Domain configuration in WPS 


Once the domain(s) is (are) configured in namelist . wps 
the geogrid . exe can be executed in order to process 
static data. 


Typing: 

. / geogrid . exe 

should start geoegrid and generate a set of 
geo em . dOX . nc files in the number corresponding to the 
number of domains configured in namelist .wps 
(paramter max_dom) 

Open generated geo_em.dOX.nc files and verify that the data 
got processed correctly (check NFUEL CAT, ZSF etc) 



Processing meteorological data 



In real cases WRF need meteorological data in order to generate initial and 
boundary conditions: 

wr f inpu t_dO 1 
wrf input d02 


initial conditions 


wrf inpu t__dOX 

wrfbdy_d02 boundary conditions 


Preparation of the meteorological data is done in two steps: 

1. Ungribing meteorological files provided in grib or grib2 format 
which is done by ungrib . exe 

2. Horizontal interpolation of the meteorological data to WRF 
domains defined by geo_em . dOX . nc files generated 
earlier which is performed by metgrid . exe . 


*m 

Processing meteorological data 

Before the ungrib can decompress grib meteorological files they must 
be linked to the location where ungrib . exe is running. 

The linking is performed using the linkgrib . sh script: 

. /link_grib . csh path_to_your_grib_f iles 

This command should generate links to your grib files in a form: 

GRIBFILE.AAA -> your_metf ilel . grb2 
GRIBFILE.AAB -> your_metf ile2 . grb2 
GRIBFILE.AAC -> your_metf ile3 . grb2 

Before running ungrib.exe make sure the correct Vtable is linked to 
your directory. For instance, if you are using NARR data you link ./ 

ungrib /Variable_Tables / Vtable . NARR toyOUrWPS directory by 
typing: 

In -s . /ungrib/Variable_Tables/Vtable . NARR ./Vtable 

you should have now in your WPS directory: 

Vtable -> . /ungrib/Variable Tables/Vtable . NARR 



*m 

Processing meteorological data 

Once the met files and Vtable are linked you can ready to 
ungrib your met data (create intermediate met files in WPS 
format). The name of the output intermediate files are defined 

in the ungrib section of namelist . wps 

&ungrib 

out_format = 1 WPS 1 , 
prefix = ' NARR_b ' , 

! prefix = ' NARR_a', 

/ 

Executing: . /ungrib . exe will start ungribing process which will generate 
intermediate met files: 

NARR_a :2012-06-2 3_1 8 
NARR_a: 2012-06-24_00 
NARR a : 2012-06-26 06 



*m 

Processing meteorological data 

Once the intermediate met files are generated, they are ready to 
be processed further by metgrid . exe, 

Executing: 

. /metgrid. exe 

should generate a list of met em.dOx files: 

met_em. dOl . 2012-06-23_00 : 00 : 00 . nc 
met_em. dOl . 2012-06-23_06 : 00 : 00 . nc 
met em. dOl . 2012-06-23 12:00:00. nc 



The next step is generating the WRF input files from met_em files using 

real.exe in . /wrf-f ire/WRFV3/test/em real/ directory 


Processing meteorological data 

Before real . exe is executed, all simulation parameters must 
be set in namelist . input, similarly as it was done for WPS 

in namelist . wps 

Once the simulation parameters are set, the WRF input files are generated 
frommet_emfilesusingreal.exe in . /wrf-f ire/WRFV3/test/ 
em_real/ directory. 

After execution of real . exe the following set of files should be generated: 

wrf input_d01 
wrf input_d02 

wrf input d03 the next step after generating these files is configuring 
— fire parameters in namelist . fire and running 

wrf. exe 



wrfbdy dOl 




/ 1 <z)i 


Configuration of fire paTa meters in 

namelist . fire 




&fuel_scalars 
cmbcnst = 17.433e+06, 
hfgl = 17. e4 , 

fuelmc_g = 0.18, 

!jc fuelmc_g = 0.09, 
fuelmc_c = 1.00, 
nfuelcats = 13, 
no_fuel_cat = 14 
/ ~ 


! scalar fuel constants 
! J/kg combustion heat dry fuel 
! W/m A 2 heat flux to ignite canopy 
! ground fuel moisture, set = 0 for dry 

! ground fuel moisture, set = 0 for dry 
! canopy fuel moisture, set = 0 for dry 
! number of fuel categories used 
! extra category for no fuel 


&fuel_categories 
fuel_name = 

'1: Short grass (1 ft) ', 

'2: Timber (grass and understory) ', 
'3: Tall grass (2.5 ft) ', 

'4: Chaparral (6 ft) ', 

'5: Brush (2 ft) ', 

'6: Dormant brush, hardwood slash', 
'7: Southern rough', 

'8: Closed timber litter', 

'9: Hardwood litter', 

'10: Timber (litter + understory) ', 
'll: Light logging slash', 

'12: Medium logging slash', 

'13: Heavy logging slash', 

'14: no fuel ' 



Configuration of tire parameters in 

namelist . fire 




fmc 

gwOl 

= 1.00000, 

0.15385, 

1.00000, 

0.31253, 

0.28571, 

0.25000, 

0.23203, 



0.30000, 

0.06625, 

0.25042, 

0.13021, 

0.11600, 

0.12065, 

0.00000, 

fmc 

gw02 

= 0.00000, 

0.07692, 

0.00000, 

0.25016, 

0.14286, 

0.41667, 

0.38398, 



0.20000, 

0.93034, 

0.16639, 

0.39149, 

0.40584, 

0.39656, 

0.00000, 

fmc 

gw03 

= 0.00000, 

0.38462, 

0.00000, 

0 . 12477, 

0.00000, 

0.33333, 

0.30801, 



0.50000, 

0.00341, 

0.41680, 

0.47830, 

0.47816, 

0.48279, 

0 . 00000, 

fmc 

gw04 

=0.0, 0.0, 

0.0, 0. 

0, 0.0, 0 

.0, 0.0, 






0.0, 0.0, 

0.0, 0 . 

0, 0.0, 0 

. 0, 0.0, 




fmc 

gw05 

= 0.00000, 

0.38462, 

0.00000, 

0.31254, 

0.57143, 

0.00000, 

0.07598, 



0.00000, 

0.00000, 

0.16639, 

0.00000, 

0.00000, 

0.00000, 

0.00000, 


/ 

Smoisture 

moisture_class_name = 
' 1 hour fuel ' , 

' 10 hour fuel ' , 

'100 hour fuel ' , 

'1000 hour fuel', 

' live fuel ' 


moisture classes= 

5, 








drying model= 

1, 

i. 

1, 

1, 

1, ! 

number of model 

- only 1= 

equilibr. 

moisture Van Wagner (1972) per 

Viney 

(1991) 

allowed 




drying lag= 

1, 

10, 

100, 

1000, 

Ie9, 

! so-called lOhr 

and lOOhr 

fuel 

wetting model= 

1, 

1, 

1, 

1, 

1, ! 




wetting lag= 

\ — i 

le9, 

le9, 

le9. 

Ie9, ! 

14 is callibrated to 


VanWagner & Pickett 1 985 , 

Canadian fire 

danger 

rating 

system 



saturation moisture^ 

2.5, 

2.5, 

2.5, 

2.5, 

2.5, ! 




saturation rain = 

8.0, 

8.0, 

8.0, 

8.0, 

8.0, ! 

(mm/h) 



rain threshold = 

LO 

O 

o 

LO 

O 

o 

LO 

O 

o 

LO 

O 

o 

0.05, ! 

mm/h rain too weak to wet 

anything 

fmc gc initialization = 

0, 

0, 

0, 

0, 

0, ! 

0 = from input. 

1 = from 

fuelmc g . 

namelist . input 2 = from 

) 

equibrium 
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Configuration of emission parameters 

in namelist . fire emissions 



Depending on the level of chemistry used appropriate namelist.fire_emissions file 
must be used. For no-chemistry runs, only passive tracer option is available. s 


To use passive tracers for smoke visualization (no chemistry) make sure you link 
appropriate emission configuration file: 

namelist . fire emissions -> namelist . fire emissions . tracers 


Remissions 


! Tracer emissions table 

! Classification : Grasslands , Grasslands , Grasslands , Open Shrublands , Open Shrublands , Open 
Shrublands , Open Shrublands , Mixed Forests , Mixed Forests , Mixed Forests , Evergreen Forest , Evergreen 


Forest , Evergreen Forest, 

compatible_chem__opt=0 , 

printsums=l, ! print sums of fuel burned and total emitted into the atmosphere 


!Fuel Cat, 

! tracer em 
trl7_l 
trl7_2 
trl7_3 
trl7_4 
trl7_5 
trl7_6 
trl7_7 
trl7 8 


_L 


LSSl 

=59, 

=59, 

=59, 

=59, 

=59, 

=59, 

=59, 

=59, 


2 , 

ons 

59, 

59, 

59, 

59, 

59, 

59, 

59, 

59, 


3 , 4 , 

[g/kg 
59, 68, 
59, 68, 
59, 68, 
59, 68, 
59, 68, 
59, 68, 
59, 68, 
59, 68, 


5 , 6 , 

burnt 
68 , 68 , 
68 , 68 , 
68 , 68 , 
68 , 68 , 
68 , 68 , 
68 , 68 , 
68 , 68 , 
68 , 68 , 


7 , 8 , 9 , 10 , 11 , 12 , 13 , 

] 

68,102,102,102,118,118,118 
68,102,102,102,118,118,118 
68,102,102,102,118,118,118 
68,102,102,102,118,118,118 
68,102,102,102,118,118,118 
68,102,102,102,118,118,118 
68,102,102,102,118,118,11 
68,102,102,102,118,118,11 


T 


i 


That is the place to define your tracer emission tactors 


Running WRF-Sfire in real cases 


• If you have: 

- generated wrf input_dOX and wrfbdy_d01 files 

- make sure the fire data are there (NFUEL_CAT, ZSF) 

- configure your meteorological and fire parameters in 

namelist . input and nemelist . fire 

- You are ready to start your simulation executing: 

. / wrf . exe 



Experiments 



1. Go to WPS dierectory 

2. Copy files to WPS folder 

3. Untar NARR, fuel and topo data 

4. Run geogrid.exe (generation of geo_em.d01.nc ...) 

5. LinkVtable 

6. Link NARR_a met data 

7. Link NARR_b met data 

8. Run ungrib.exe (generation of NARR_a/b interm, files) 

9. Run metgrid.exe (generation of met_em.d01_xx files) 

10. Run real.exe (generation of wrfinput and wrfbdy files) 

11. Run wrf.exe 



Experiments 



1. Go to /scratch/wrf-f ire/WPS 

2. Copy files to WPS folder: 

cp /home/netapp-clima-scratch/smr2 7 17 /lab_sess ions /Experiment 6 /namelist .wps . 
cp /home/netapp-clima-scratch/smr2717/lab_sessions/Experiment6/namelist . input . 
cp /home/netapp-clima-scratch/smr2717/lab_sessions/Experiment6/NARR. tgz . 
cp /home/netapp-clima-scratch/smr2717/lab_sessions/Experiment6/landf ire_data. tgz . 
cp /home/netapp-clima-scratch/smr2717/lab_sessions/Experiment6/ned_data. tgz . 
cp /home/netapp-clima-scratch/smr2717/lab_sessions/Experiment6 /GEOGRID . TBL . 

3. Untar NARR fuel and topo data 

tar — xvf NARR. tgz 

tar —xvf landf ire_data. tgz 

tar —xvf ned_data.tgz 

4. Run . /geogrid. exe 
it should generate: 

geo_em. dOl . nc geo_em. d02 . nc geo_em. d03 . nc geo_d04.nc 

5. Link Vtable: 


In -s . /ungrib/Variable_Tables/Vtable .NARR ./Vtable 

6. Link NARR_a met data 

./linkgrib ./NARR_a/* 



Experiment6 

6. Link NARR_a met data 

./linkgrib . /NARR_a/* 



That will generate a lot of GRIBFILE links: 


GRIBFILE.AAA GRIBFILE. AAE 
GRIBFILE . AAB GRIBFILE . AAF 
GRIBFILE . AAC GRIBFILE . AAG 
GRIBFILE . AAD GRIBFILE . AAH 


GRIBFILE. AAI 
GRIBFILE. AAJ 
GRIBFILE. AAK 
GRIBFILE. AAL 


GRIBFILE. AAM 
GRIBFILE. AAN 
GRIBFILE. AAO 
GRIBFILE. AAP 


7. Modify namelist.wps (f.e emacs namelist . wps) to have: 

end_date = ' 2012-06-24_18 : 00 : 00 ' , ' 2012-06-24_18 : 00 : 00 ' , 

' 2 012-0 6-2 4_1 8 : 00 : 00 ' , ' 2012-06-24_18 : 00 : 00 ' , 

8. Run ungrib to extract the met data 
. /ungrib. exe 

That will generate ungribed NARR_a files: 

NARR_a :2012-06-2 3_1 8 
NARR_a :2012-06-2 4_0 0 
NARR a : 2012-06-26 06 



Experiment6 



7. Remove old GRIBFILE links: rm GRIBFILE.* and link NARR b met data 


. /linkgrib . /NARR_b/ * 

That will generate again a lot of GRIBFILE links linked to NARR_b: 


GRIBFILE. AAA 
GRIBFILE. AAB 
GRIBFILE. AAC 
GRIBFILE. AAD 


GRIBFILE. AAE 
GRIBFILE. AAF 
GRIBFILE. AAG 
GRIBFILE. AAH 


GRIBFILE. AAI 
GRIBFILE. AAJ 
GRIBFILE. AAK 
GRIBFILE. AAL 


8. Edit namelist. wps 


emacs namelist. wps 


change : 


! prefix = ' NARR_b ' , 
prefix = ' NARR_a ' , 
to 


GRIBFILE. AAM 
GRIBFILE. AAN 
GRIBFILE. AAO 
GRIBFILE. AAP 


prefix = ' NARR_b ' , 

! prefix = 'NARR_a', 

9. Run ungrib again to extract the second set of meteo data 
. /ungrib. exe 

That will generate ungribed NARR_b files: 

NARR_b :2012-06-2 3_1 8 
NARR_b: 2012-06-24_00 
NARR b : 2012-06-26 06 



Experiments 



10. Link METGRID.TBL to WPS directory: 

In -s . /metgrid/METGRID . TBL ./ 

11. Run metgrid.exe 
. /metgrid.exe 

That will generate a lot of met_em.dOX.xxxx files 

12. Goto /wrf-f ire/WRFV3/test/em_real/ , and copy namelist . input there 
cd /wrf-f ire/WRFV3 /test /em_real/ 

cp ../../.. /WPS/namelist . input ./ 

13. link met_em files to /wrf-f ire/WRFV3/test/em_real/ 

In -s ../../.. /WPS/met_em. * ./ 

14. Edit namelist. input to set the simulation end time to: 


end_year 

= 2012, 

2012, 

2012, 

2012 

end_month 

= 06, 

06, 

06, 

06, 

end_day 

= 24, 

24, 

24, 

24, 

end hour 

= 06, 

06, 

06, 

06, 


15. Run ./real . exe to generate 

wrf input_d01 , wrf input_d02 , wrf inout_d03 , wrf input_d04 , 
wrfbdy_d01 

16. Run wrf . exe 


